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Abstract 

Positive frequency Wightman function, vacuum expectation values of the field square 
and the energy-momentum tensor induced by a Z2-symmetric brane with finite thickness 
located on {D + 1)- dimensional AdS background are evaluated for a massive scalar field 
with general curvature coupling parameter. For the general case of static plane symmetric 
interior structure the expectation values in the region outside the brane are presented as 
the sum of free AdS and brane induced parts. For a conformally coupled massless scalar 
the brane induced part in the vacuum energy-momentum tensor vanishes. In the limit of 
strong gravitational fields the brane induced parts are exponentially suppressed for points 
not too close to the brane boundary. As an application of general results a special model 
is considered in which the geometry inside the brane is a slice of the Minkowski spacetime 
orbifolded along the direction perpendicular to the brane. For this model the Wightman 
function, vacuum expectation values of the field square and the energy-momentum tensor 
inside the brane are evaluated as well and their behavior is discussed in various asymptotic 
regions of the parameters. It is shown that for both minimally and conformally coupled 
scalar fields the interior vacuum forces acting on the brane boundaries tend to decrease the 
brane thickness. 



PACS: 0A.62.+V; ll.lO.Kk 



1 Introduction 

Recent proposals of large extra dimensions use the concept of brane as a sub-manifold embedded 
in a higher dimensional spacetime, on which the Standard Model particles are confined (see, for 
instance [Ij). Braneworlds naturally appear in string/M-theory context and provide a novel 
setting for discussing phenomenological and cosmological issues related to extra dimensions. 
The models introduced by Randall and Sundrum are particularly attractive. In RS 2-brane 
model [2j the corresponding background solution consists of two parallel flat 3-branes, one with 
positive tension and another with negative tension embedded in a five dimensional AdS bulk. 
The fifth coordinate is compactified on /Z2, and the branes are on the two fixed points. It is 
assumed that all matter fields are confined on the branes and only the gravity propagates freely 
in the five dimensional bulk. In RS 1-brane model [3j there is only one positive tension brane in 
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a five dimensional AdS bulk. This model is obtained from 2-brane one by sending the negative 
tension brane to infinity. Due to the warp factor in the AdS metric the infinite extra dimension 
gives a finite contribution to the volume. More recently, alternatives to confining particles on 
the brane have been investigated and scenarios with additional bulk fields have been considered. 

Motivated by the problems of the radion stabilization and the generation of cosmological 
constant, the role of quantum effects in braneworlds has attracted great deal of attention [l]-|42j. 
A class of higher dimensional models with compact internal spaces is considered in |43j . Many 
of treatments of quantum fields in braneworlds deal mainly with the case of the idealized brane 
with zero thickness. This simplification suffers from the disatvantage that the curvature tensor is 
singular at the brane location. In addition, the vacuum expectation values of the local physical 
observables diverge on the brane. From a more realistic point of view we expect that the branes 
have a finite thickness and the thickness can act as natural regulator for surface divergences. 
The finite core effects also lead to the modification of the Priedmann equation describing the 
cosmological evolution inside the brane. In particular, a phase with accelerated expansion on 
the brane can be obtained in the model with pressureless matter as the only non-gravitational 
source in the modified Friedmann equation. In sting theory there exists the minimum length 
scale and we cannot neglect the thickness of the corresponding branes at the string scale. The 
branes modelled by field theoretical domain walls have a characteristic thickness determined by 
the energy scale where the symmetry of the system is spontaneously broken. Various models 
are considered for a thick brane (see, for instance, [44] and references therein). Mainly, these 
models are constructed as solutions to the coupled Einstein-scalar equations by choosing a 
suitable potential for the scalar field. Vacuum fiuctuations for a simple thick de Sitter brane 
supported by a bulk scalar field with an axion like potential and the self-consistency of this 
braneworld are investigated in Refs. |40j . 

In the present paper we investigate the effects of core on properties of the quantum vacuum for 
the general plane symmetric static model of the brane with finite thickness. The most important 
quantities characterizing these properties are the vacuum expectation values of the field square 
and the energy-momentum tensor. Though the corresponding operators are local, due to the 
global nature of the vacuum, the vacuum expectation values describe the global properties of 
the bulk and carry an important information about the internal structure of the brane. In 
addition to describing the physical structure of the quantum field at a given point, the energy- 
momentum tensor acts as the source of gravity in the Einstein equations. It therefore plays an 
important role in modelling a self-consistent dynamics involving the gravitational field. As the 
first step for the investigation of vacuum densities we evaluate the positive frequency Wightman 
function for a massive scalar field with general curvature coupling parameter. This function 
gives comprehensive insight into vacuum fiuctuations and determines the response of a particle 
detector of the Unruh-DeWitt type moving in the brane bulk. The problem under consideration 
is also of separate interest as an example with gravitational and boundary-induced polarizations 
of the vacuum, where all calculations can be performed in a closed form. The corresponding 
results specify the conditions under which we can ignore the details of the interior structure and 
approximate the effect of the brane by the idealized model. In addition, as it will be shown 
below, the phenomenological parameters in the zero-thickness brane models such as brane mass 
terms for scalar fields are calculable in terms of the inner structure of the brane within the 
framework of the model considered in the present paper. 

The plan of this paper is as follows. In Section [2] we consider the Wightman function in 
the exterior of the brane for the general structure of the core assuming that the components 
of the metric tensor and their derivatives are continuous at the transition surface between the 
core and the exterior. The Section [3] is devoted to the generalization of the corresponding 
results when an additional surface shell is present on the bounding surfaces between the core 
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and the exterior. By using the formulae for the Wightman function, in Section 2] we investigate 
the vacuum expectation values of the field square and the energy-momentum tensor. As an 
illustration of the general results, in Section [S] we consider a model with Minkowskian geometry 
inside the brane. For this model the vacuum expectation values inside the core are investigated 
as well. The last section contains a summary of the work. 

2 Wightman function 
2.1 Model 

We consider a domain wall with finite thickness on background of {D + l)-dimensional AdS 
spacetime. As in RS 1-brane scenario we will assume that our model is Z2-symmetric with 
respect to the plane y = located at the center of the brane. The spacetime is described by 
two distinct metric tensors in the regions outside and inside the brane. The corresponding line 
element in the exterior region |y| > a has the form 

ds'^ = gikdx'dx'' = e'^^^^y^rj^udx^'dx'' - dy'^, (2.1) 

where r/^^ = diag(l, —1, . . . , —1) is the metric for the D-dimensional Minkowski spacetime, l/kri 
is the AdS curvature radius. Here and below i, k = 0,1, . . . , D, and //, = 0, 1, . . . , D — 1. We 
will assume that inside the brane (region \y\ < a) the spacetime geometry is regular and is 
described by the static plane symmetric line element 

ds^ = ^<y)df - e^<y'^d^^ - e^^^y^y^ (2.2) 

where x = (x^, . . . ,x^~^) are coordinates parallel to the brane. Due to the Z2-symmetry the 
functions u{y), v{y), w{y) are even functions on y. These functions are continuous at the core 
boundary: 

u{a) = v{a) = —k^a, w{a) = 0. (2-3) 

Here we assume that there is no surface energy-momentum tensor located at \y\ = a and, hence, 
the derivatives of these functions are continuous as well. The generalization to the case with an 
infinitely thin shell at the boundary of two metrics will be discussed in the next section. 

For the metric corresponding to line element ()2.2p the nonzero components of the Ricci tensor 
are given by expressions (no summation over i) 

e-2- ^u" + n'2 + {D- l)u'v' - u'w'] , 

^-2w ^^// + (2) - l)v'^ + v'u' - v'w'] , i = l,...,D-l, (2.4) 
e-S"- + {D- l)v" + n'2 + {D - l)v'^ - w'v! - {D - l)w'v'] , 

where the prime means the derivative with respect to the coordinate y. The corresponding Ricci 
scalar has the form 

R = -2e^2«- [^// + (2) _ i)v" + + D{D - l)v'^/2 

+{D - l)uv' - u'w' -{D- l)v'w'] . (2.5) 

In the special case u{y) = v{y) one has Rq = R\, i = 1, . . . ,D — 1, and the model is Poincare 
invariant along the directions parallel to the brane. 
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2.2 Eigenfunctions 

In this paper we are interested in the vacuum polarization effects for a scalar field with general 
curvature coupling parameter ^ propagating in the bulk described above. The corresponding 
field equation has the form 

(V.V* + m^ + ^R)^ = 0, V.V* = -^di{g'''^/\i\dk), (2.6) 

where Vj is the covariant derivative operator associated with line element (j2.ip outside the brane 
and with line element (12. 2p inside the brane. The values of the curvature coupling parameter 
^ = and ^ = with = (D — 1)/4Z) correspond to the most important special cases of 
minimally and conformally coupled scalar fields. As a first stage for the evaluation of the vacuum 
expectation values (VEVs) for the field square and the energy-momentum tensor we consider 
the positive frequency Wightman function {0\ip{x)ip{x')\0) , where |0) is the amplitude for the 
corresponding vacuum state. This function also determines the response of the Unruh-DeWitt 
type particle detector at a given state of motion (see, for instance, [l5]). The Wightman function 
can be evaluated by using the mode sum formula 

{0\^ixMx')\0) =^^^ix)ipUx'), (2.7) 

a 

where {{pa{x),(p'^{x')} is a complete orthonormalized set of positive and negative frequency 
solutions to the field equation. The collective index a can contain both discrete and continuous 
components. In Eq. (j2.7p it is assumed summation over discrete indices and integration over 
continuous ones. 

On the base of the plane symmetry of the problem under consideration the corresponding 
eigenfunctions can be presented in the form 

^aix') = M^nixiy), (2.8) 

where the standard Minkowskian modes on the hyperplane parallel to the brane: 



</.k(x^) = ^====, u; = VWT)^, k = \k\, (2.9) 

A is the separation constant. Below we will assume that y ^ 0. The corresponding formulae in 
the region y < are obtained from the Z2-symmetry of the model. Substituting eigenfunctions 
(12.81) into field equation (12. Op . for the function f{y) one obtains the following equation 

^-u~iD-i)v-wQ^ Lu+iD-i)v-wQj] _ ^^2 + + k^e-^" - e-2") - A^e-^^] fx = 0. (2.10) 

For the exterior AdS geometry one has u{y) = v{y) = —koy, and R = —D{D + l)k'jj. In this 
case the solution to equation (|2.10p for the region y > a is 

fx{y) = e'^'^y/^ [AMXz] + BMXz)] , z = e^^^kn, (2.11) 

where Ai, and Bi, are integration constants, Juix), Yiy{x) are the Bessel and Neumann functions, 
and we use the notation 

D'^/A- D{D + 771"^/ kj). (2.12) 

Here we will assume values of the curvature coupling parameter for which is real. For imaginary 
v the ground state becomes unstable 06j. Note that for a conformally coupled massless scalar one 
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has V = \/2 and the cyhnder functions in Eq. (|2.11|) are expressed via the elementary functions. 
The regular solution of equation (|2.10p in the region y < a we will denote by R{y,k, X). Note 
that the parameters k and A enter in the radial equation in the form fe^ and A^. As a result 
the regular solution can be chosen in such a way that R{y, —k, A) = const • R{y, k, A) and 
R{y, k, — A) = const • R{y, k, A). Now for the radial part of the eigenfunctions one has 



f(\-) ^(y'^>-^)> fory<a , , 

h[y) - i ^Dk,y/2 [AM^z) + i?.y.(Az)] , for y > a ■ ^^"^^^ 



The coefficients Ai, and Bi, are determined by the conditions of continuity of the radial function 
and its derivative at y = a. Prom these conditions we find 

= e'"'^''/^ [AMXza) + BMXza)] , (2.14) 
K = kDe'"'^''^^A,[DMXza)/2 + XzaJl{Xza)] 

+B, [DY,(Xza)l2 + XzaYXXza)] ] , (2.15) 

where 

Ra = R{a, k, A), < = dyR{y, k, X)\y=a. (2.16) 

Solving these equations with respect to A^ and B^, and using the Wronskian for the Bessel and 
Neumann functions, one has 

A, = |e-^^°"/2i?,y,(AzJ, B, = -le-'^'^'^/^RaMXza). (2.17) 

Here and in what follows for a cylinder function F{z) we use the notation 

m . + (f - ^) F(.). (2.18) 

Note that due to our choice of the function R{y,k,X), the logarithmic derivative in formula 
(I2.18P is an even function on A. Hence, in the region y > a the radial part of the eigenfunctions 
has the form 

fx{y) = \e''^^^y-'^'^'^Ra9u{Xza, Xz), (2.19) 

where the notation 

g,{u, v) = Y,{u)Mv) - Mu)Y,{v) (2.20) 

is introduced. 

For the eigenfunctions we have the following orthonormalization condition 



dyd^-'xy'\g\g^\^{x)vl,{x) = -j^, (2.21) 

where S^a' is understood as the Kronecker symbol for discrete indices and as the Dirac delta 
function for continuous ones. Substituting eigenfunctions ()2.8p . the normalization condition is 
written in terms of the radial eigenfunctions 

+00 

dy^A9\9'''fx{y)fx'iy) = Six - A'). (2.22) 

■oo 

As the integral on the left is divergent for A' = A, the main contribution in the coincidence 
limit comes from large values \y\. By using the expression ()2.19p for the radial part in the 
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region \y\ > a and replacing the Bessel and Neumann functions by the leading terms of their 
asymptotic expansions for large values of the argument, we obtain the following result: 



This relation determines the normalization coefficient for the interior eigenfunctions. 

In addition to the modes with real A discussed above, modes with purely imaginary A can 
exist. For this type of modes the radial eigenfunctions in the region y > a are given by the 
expression C^e^^^^^'^ Ku{\X\z) , where Ky{x) is the MacDonald function. From the continuity of 
the eigenfunctions at y = a one finds Ra = C^e^^^"/^i^jy(| Ajza), and from the continuity of the 
radial derivative we see that |A| is a solution of the equation 

K{\\\za) , D _ dyR{y,k,e^'lM)\y=a 



£(|A|z) + 2 - ' kDR{a,k,e-^^/^\\\) ' ^^'^^^ 

As it follows from (|2.9p . for the modes with purely imaginary A the corresponding frequency 
is imaginary in the region < |A| and the vacuum state becomes unstable. To have a stable 
ground state, in the discussion below we assume that there are no modes with imaginary A. 

2.3 Wightman function in the exterior region 

Substituting the eigenfunctions (|2.8|) into the mode sum (|2.7p . for the Wightman function one 
finds 

uD-l p ik-(x-x') poo \ 

(0|^(x)v.(x')|0) = ^{zz')- / dk— _ / dA- 



j2(Az„)+y,2(AzJ 
For further transformation of this formula we use the identity 



(2.25) 



g.(Az A.MA. AzO ^ j^^^^^j^^^,^ _ 1 ^ J^f^yj)^Xz)Hi^){Xz\ (2.26) 

where Hy^\x), s = 1,2 are the Hankel functions. This allows us to present the Wightman 
function in the form 

(0|(^(x)(^(x')|0) = ^(05|v^(x)<^(xO|05> + ((^(x)(^(x'))b, (2.27) 



where 



uD-l p ik-(x-x') poo \ 

(05|^(x)(^(x')|05) = -^{zz')- / dk-— ^ / dX- 



2 ' ' J (27r)^-i io ^/WTX? 
xJ^{Xz)MXz')e''^''-'^^''"+^\ (2.28) 

is the positive frequency Wightman function for the AdS spacetime without boundaries (see, for 
instance, [Ml)) and the part 

l-D-l p ik-(x-x') 2 „oo 

{ifixMx'))^ = -^(zz')- J dk^^-^Y^J^ dXX 

X y-;""\ ffW(A.)ffW(A.O, (2.29) 

VA2 + k^ H)f\XZa) 
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induced by the brane. Quantum effects in free AdS spacetime are well investigated in literature 
(see references given in |34] ) and in the discussion below we will be mainly concentrated on the 
effects induced by the brane. 

On the right of formula (j2.29p we rotate the integration contour in the complex plane A by the 
angle -7r/2 for s = 1 and by the angle —tt/2 for s = 2. Under the condition z + z' > 2za + \t — t'\, 
the contributions from the integrals over the arcs of the semicircle with large radius in the 
upper/lower half-plane tend to zero. In particular, this condition is satisfied in the coincidence 
limit for points outside the brane. Further, by using the property that the logarithmic derivative 
of the function R{y,k, X) in formula (I2.18P is an even function on A, we can see that the integrals 
over the segments (0, ik) and (0, —ik) of the imaginary axis cancel out. As a result, after 
introducing the modified Bessel functions, the brane- induced part can be presented in the form 

ivixMx'))^ = -^^{zz')T [dk(^'^<^-^'^ TdW^"^^'''^ 



{27:)D'-' J — """'K,{XZa) 
K^{\z)K^{\z') 



VA2 - fc2 



■ cosh 



y/X^ -k'^{t' -t) . (2.30) 



Here and below the tilted notation for the modified Bessel functions Iu{x) and Ky{x) is defined 
by the formula 

F{x) = xF'{x) + TZ{a, k, x)F{x), (2.31) 

with the notation 

, , , L» dyR(y,k,xe'^'/'^/za)L-a 
7^ a, k,x) = -- \ ' W2/ V • 2.32 
2 k£)R[a, k,xe^^'^ / Za) 

Note that, by taking into account this notation, equation (|2.24p for the modes with imaginary 
A is written in the form K,y{\X\za) = 0. 

As we see from (j2.30p , the information about the inner structure of the brane is contained 
in the logarithmic derivative of the interior radial function in formula (j2.32p . In formula (|2.30p 
the integration over the angular part can be done by using the formula 

(27r)(^-i)/2 " io ^ ^ (fc|x|)(^-3)/2 ' i^-^-^^ 

for a given function F{k). In the RS 1-brane model with the brane of zero thickness the 
brane induced part in the Wightman function is given by a formula similar to (j2.3Up with the 
replacement [33] 

n{a,k,x) ^ D/2-2DC- c/2kD, (2.34) 

in the definition (j2.3ip of the tilted notation. On the right of relation (j2.34p . the parameter c 
is the brane mass term for a scalar field which is a phenomenological parameter in the model 
with zero thickness brane. As we see, in the model under consideration the effective brane mass 
term is determined by the inner structure of the core. Note that in RS 2-brane model the mass 
terms on separate branes determine the 1-loop effective potential for the radion field and play 
an important role in the stabilization of the interbrane distance. 



3 Model with an infinitely thin shell on the brane boundary 

The results considered in the previous section can be generalized to the models where an addi- 
tional infinitely thin plane shell located at y = a is present with the surface energy-momentum 
tensor rf^. Let n* be the normal to the shell normalized by the condition n^n* = — 1. We assume 
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that it points into the bulk on both sides of the shell. From the Israel matching conditions one 
has 

{Kik-Khik} = 87rGT,k, (3.1) 

where the curly brackets denote summation over each side of the shell, /ij^ = gik + ^jn^ is the 
induced metric on the shell, Kik = h'^h^^Vrng its extrinsic curvature and K = K^. For the shell 
located on the boundary y = a and for the region y ^ a one has = 5^e"''^^-* and the non-zero 
components of the extrinsic curvature are given by the formulae 

= -n'(a-), Ki = -djv'ia-), i = 1,2, , . . . , D - 1. (3.2) 

The corresponding expressions for the region y ^ a are obtained by taking u{y) = v{y) = —kay, 
w{y) = and changing the signs for the components of the extrinsic curvature tensor: 

Kl = -6{kD, i = l,2,,...,D-l, y = a + . (3.3) 

Now from the matching conditions (13. Ih we find (no summation over i) 

u'{a-) = -ko + SttG (^i - , i = l,2,...,D-l, (3.4) 

v'{a-) = -ko + ^^rl (3.5) 

where f'{a—) is understood in the sense limy_+a-o /'(y)- For models with Poincare invariance 
along the directions parallel to the brane one has = Tq . 

The discontinuity of the functions u'{y) and v'{y) at y = a leads to the delta function term 

2 [u'ia-) + {D- l)v'{a-) + Dko] S{y - a) (3.6) 

in the Ricci scalar and, hence, in the equation (j2.10p for the radial eigenfunctions. Note that 
the expression in the square brackets is related to the surface energy-momentum tensor by the 
formula 

u'{a-) + {D-l)v'{a-) + DkD = ^^T, (3.7) 

where r is the trace of the surface energy-momentum tensor. 

For a non-minimally coupled scalar field, due to the delta function term in the equation 
for the radial eigenfunctions, these functions have a discontinuity in their slope at y = o. The 
corresponding jump condition is obtained by integrating the equation (|2.10p through the point 
y = a: 

/;(«+) - f'xia-) = ^-^rhia). (3.8) 

Now the coefficients in the formulae ()2.13p for the eigenfunctions are determined by the continu- 
ity condition for the radial eigenfunctions and by the jump condition for their radial derivative. 
It can be seen that the corresponding eigenfunctions are given by the same formulae (j2.19p and 
(I2.23P with the new barred notation 

F(.)..f'W+(--^^^r-^)fW. (3.9) 

Hence, the part in the exterior Wightman function induced by the brane is given by formula 
()2.30p . where the tilted notation is defined by Eq. (j2.3ip with the function 

, , D levrG^ dyR(y,k,xe'''/'^/za)L=a 
n{a, k,x) = -- I T - \ ' J "/'V - 3-10 

2 (iJ — kDR{a,k,xe^^i^ / Za) 
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The trace of the surface energy-momentum tensor in this expression is related to the components 
of the metric tensor inside the brane by formula (j3.7p . As before, in the model with surface 
energy-momentum tensor the equation for the modes with imaginary A is written in the form 

k,{\\\Za)=Q. 



4 Vacuum expectation values outside the brane 
4.1 VEV of the field square 

The VEV of the field square is obtained by computing the Wightman function in the coincidence 
limit x' — > X. In this limit expression (j2.27p gives a divergent result and some renormalization 
procedure is needed. Outside the brane the local geometry is the same as that for the AdS 
spacetime. Hence, in the region y > a the renormalization procedure for the local characteristics 
of the vacuum, such as the field square and the energy-momentum tensor, is the same as for the 
free AdS spacetime. 

By using the formula for the Wightman function from the previous section, the VEV of the 
field square in the exterior region is presented in the form 

{0\^'\0) = l{0s\^'\0s) + {v\, (4.1) 

where {Os\^'^\Os) is the VEV of the field square in the free AdS spacetime. The part induced 
by the brane is directly obtained from formula (|2.30p in the coincidence limit: 

2^-ir((Z) - l)/2) io A ^^^i^.(Az.)^/A^• ^'"'^ 

The VEV of the field square in the free AdS spacetime is well investigated in literature [17] and 
is given by the formula 

This VEV does not depend on the spacetime point, which is a direct consequence of the maximal 
symmetry of the AdS bulk. For even D, the expression on the right of (j4.3p is finite and can be 
taken as the renormalized value for the field square. For odd D, an additional renormalization 
removing the pole of the gamma function in the numerator is necessary. 

By taking into account that the tilted functions in formula (j4.2p depend on A;, in the general 
case of the core model the corresponding expression cannot be further simplified. This can be 
done in the model with Poincare invariance along the directions parallel to the brane. For this 
special case in (j2.2p we have u{y) = v{y) and the term with k'^ in (j2.10p disappears. As a result 
the interior radial functions R{y, k, A) do not depend on k, R{y, k, A) = R{y, A), and, hence the 
same is the case for the tilted functions in ()4.2p . In this case the formula for the VEV of the 
field square is further simplified by using the formula 

' .tk"-^ r f^L = ^ix' /" (4.4) 



'0 Jk 

This leads to the result 



VA2 - fe2 2r (f ) 



= -T^^i^^^ r dxx^-' ¥^K'Axz), (4.5) 



(47r)^/2r(Z?/2)7o K^xza) 
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where the tilted notations are defined by (|2.31|) . 

The integral in Eq. (j4.5p is exponentially convergent in the upper limit for z > Za and 
diverges for points on the boundary of the brane. The surface divergences in the VEVs of local 
physical observables bilinear in the field are well known in quantum field theory with boundaries. 
They can be regularized by considering more realistic model with smooth transition between the 
interior and exterior metrics. By taking into account that for a given z the main contribution 
into the integral in (j4.5p comes from the region x < — Za)~^, we expect that under the 
condition y — a » yo, with yo being the thickness of the transition range, the results of the 
present paper will be valid for this kind of models as well. 

At large distances from the brane, z ^ Za, we introduce a new integration variable y = xz 
and expand the integrand over Za/z. By using the formula for the integral involving the square 
of the MacDonald function, to the leading order we obtain 

2. k^~\zjzf'^ 7^(a) + z.^(Zj/2 + z.)^(Z)/2 + 2^) 

W ;b 2^+2'^+i7r(«-i)/2 7^(a) - u i'r^{u)T{{D + l)/2 + z^) ' ^ ^ 

where we have introduced the notation 

'^^''> ~ 2 {D-l)kD kDR{a,0) ■ ^ 

As we see, at large distances from the brane the brane induced part is exponentially suppressed 
by the factor exp{—2ukDy)- 

4.2 Vacuum energy-momentum tensor 

Now we turn to the investigation of the VEV of the energy-momentum tensor in the region 
y > a. Having the Wightman function and the VEV for the field square, this can be done on 
the base of the formula 



{0\Tik\0) = lim did'f,{0\^ixMx')\0) + 



(0|^2|o). (4.8) 



Note that on the left of this formula we have used the expression for the energy-momentum 
tensor of a scalar field which differs from the standard one by the term which vanishes on the 
solutions of the field equation (12. 6p (see Ref. [48]). Similar to the Wightman function, the 
components of the vacuum energy-momentum tensor are presented in the decomposed form 

{OlTiklO) = l{0s\Tik\0s) + {Tik)h, (4.9) 

where {Os\Tik\Os) is the vacuum energy-momentum tensor in the free AdS spacetime and the 
part {Tik)h is induced by the brane. As the VEV of the field square in free AdS spacetime does 
not depend on spacetime point, the corresponding VEV is found from the formula {Os\Tj^\Os) = 
6^771^ {Osl^'^l^s) by taking into account formula (j4.3p . For a conformally coupled massless scalar 
field and for even values D the renormalized free AdS part in the VEV of the energy-momentum 
tensor vanishes. For odd values of D, this part is completely determined by the trace anomaly 
(see, for instance, [45]). 

Substituting the expressions of the Wightman function and the VEV of the field square into 
formula ()4.8p . for the part of the energy-momentum tensor induced by the brane one obtains 

.k, _ ^E+'^'^^f r ^^ ^ ^-(^^'^) ^^^^ [k,KA^z)] 



{Tt)h = — TT^rirr^hf^^nr-rT dkk""'' dAA - — ' (4.10) 
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where for a given function g{v) we use the notations 

F^'Hk,9{v)] = Q-2e) [t;V'W + (i^+4^)^5WW (4.11) 

+(v' + i,^)g\v)]+iz^k^-v')g\v) 
F^)%g{v)] = F(^)[k,g{v)] + [v^-^^^^g\v), i = l,...,D-l, (4.12) 

F^^\k,g{v)] = -f.g'\v) + £(i^-l)vg(v)g'iv) 

+\ [v"^ + + 2e^(^ + 1) - ^V2] g\v). (4.13) 

As in the case of the field square, these expression may be further simphfied in the model with 
Poincare invariance along the directions parallel to the brane for which u{y) = v{y). In this case 
the tilted functions in (j4.10p do not depend on k and after the integration one finds 

where for a given function g{v) we have introduced the notations 



v^g'\v)+[D + ^:^]vg{v)g'{v) 



F^'^hgiv)] = -^g'\v) + j{AC-l)vg{v)g'iv) 



0,1,..., D-1, (4.15) 



+\ [v"^ + ^'^ + 2e^(^ + 1) - ^V2] g'^iv). (4.16) 

It can be explicitly checked that the VEVs given by ()4.14p satisfy the continuity equation 
Vfc(T^^)b = 0, which for the geometry under consideration takes the form 

z^+ia,(z-^(rii')b) + D{T^)^ = 0. (4.17) 

For a conformally coupled massless scalar field one has = 1/2 and from formulae (I4.15|) . ()4.16p 
it follows that F^^[K^{x)] = F^^\Ky{x)] = 0. Hence, in this case the brane induced parts in 
the VEVs of the energy-momentum tensor vanish. Note that for a conformally coupled scalar 
and for even values D the conformal anomaly is absent and the free AdS part in the vacuum 
energy-momentum tensor vanishes as well. 

For large distances from the brane, z ^ z^, introducing a new integration variable y = xz 
we expand the integrand over Za/z. To the leading order this leads to the result 

The integrals in this formula may be evaluated by using the formulae from ^49j . Note that the 
free AdS parts in the VEVs of both field square and the energy-momentum tensor do not depend 
on the spacetime point and, hence, at large distances from the brane they dominate in the total 
VEVs. Noting that in the limit of strong gravitational field in the region outside the brane, 
corresponding to large values ko, one has z/za = e'^^'^y-"^ > 1, we see that formulae (j4.6p . 
(|4.18p also describe the asymptotic behavior of the brane induced VEVs in this limit. Hence, 
in the limit of strong gravitational field, for the points not too close to the brane boundary, the 
brane induced parts are exponentially suppressed. The free AdS parts behave as k^~^ and their 
contribution dominates for strong gravitational fields. 



11 



5 Model with flat spacetime inside the brane 



As an application of the general results given above let us consider a simple example assuming 
that the spacetime inside the brane is flat. The corresponding models for the cosmic string and 
global monopole cores were considered in Refs. j-lOj and are known as flower-pot models. As we 
will see below, this model is exactly solvable in the sense that the VEVs inside the brane are 
expressed in terms of closed analytic formulae. This model is also of interest from a physical 
point of view as it reduces to the standard RS 1-brane model in the limit of the zero thickness 
brane and, hence, is the simplest generalization of the latter. From the continuity condition on 
the brane boundary it follows that in coordinates the interior line element has the form 

ds'^ = e-^^'^^'r^^ydx^'dx'' - dy^. (5.1) 

In the rescaled coordinates 

Xf" = e-^^'^x^, (5.2) 

this line element is written in the standard Minkowskian form. From the matching conditions 
(13. 4p . (13. 5p we find the corresponding surface energy-momentum tensor with the non-zero com- 
ponents 

= ^^fcz5<5f, i = 0,l,...,D-l. (5.3) 

The corresponding surface energy density is positive. Lets us consider the VEVs in the exterior 
and interior regions separately. 

5.1 Exterior region 

For the simple model under consideration the radial equation (j2.10p is easily solved in the interior 
region. The corresponding radial eigenfunction with Z2-symmetry, R{—y, k, A) = R{y, k, A), has 
the form 

R{y, k, A) = Cu cos{kyy), k^ = X^e^^^" - . (5.4) 
The coefficient Cy is found from formula (I2.23P 



with the barred notation for the cylindrical functions 



F{z) = zF'{z) + 



D k 



2^L> + -f!- tanikyo) 



F{z). (5.6) 



As a result, the parts in the Wightman function, in the VEVs of the field square and the 
energy-momentum tensor induced by the brane are given by formulae (|2.30p . (j4.5p and (j4.14p 
respectively, where the tilted notations for the modified Bessel functions are defined by (j2.3ip 
with the coefficient 



TZ{a,k,x) = y - 2^D - ^/x^ + m'^/kl)tanh{akD\/x'^ + m?/ky). (5.7) 



By using the recurrence relation for the function Ku{x), it can be seen that under the condition 

u + 2^D- D/2 + {m/kD)tanh{am)^0 (5.8) 

one has Ky{x) < and in the model under consideration there are no modes with imaginary A. 
In particular, this condition is satisfied for minimally and conformally coupled fields. For large 
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distances from the brane the asymptotic behavior of the parts induced by the brane is given 
by formulae (j4.6p and ()4.18p for the field square and the energy-momentum tensor respectively. 
Now in these formulae we have 

7^(a) = D/2 - 2^D - {m/ko) tanh(am). (5.9) 

Comparing (j5.7p with (|2.34p . we see that in the limt a ^ from the results of the model with 
flat interior spacetime the corresponding formulae in the RS 1-brane model with a zero thickness 
brane are obtained. 

As it was mentioned before the VEVs of the field square and the energy-momentum tensor 
diverge on the boundary of the brane. To find the leading terms in the corresponding asymptotic 
expansions for points near the boundary, we note that for these points the main contribution 
into the integrals come from large values x and we can replace the modified Bessel functions by 
corresponding asymptotic formulae for large values of the argument. In this way, assuming that 
y — a <^ a, one finds 

2 kpAn TUP - l)/2) 

with the notation 

Ai^ = AD-D^ + l + 4(D{D - 3) - Am^/kjj. (5.11) 

For D ^ 2 the VEV of the field square is finite on the core boundary. Note that in the model 
with zero thickness brane located at y = 0, the corresponding VEVs near the brane behave as 

The asymptotic behavior for the energy-momentum tensor is found in the similar way. To 
the leading order we have (no summation over i) 

knAn{^-^D)r{{D + l)/2) 
)b - - 2Z.+4^(D+i)/2 ^y_a)D ' ^ = 0,l,...,D-l. (5.12) 

The behavior for the radial stress is found most easily from the continuity equation ()4.17p . This 
leads to the result 

Dkpiy-a) q 
{-^D/h' D -I \^0/b- (5-13) 

For a conformally coupled field the most easy way to find the asymptotic behavior is to note 
that one has the following relation for the trace of the energy-momentum tensor 

From here and the continuity equation it follows that for a conformally coupled scalar field we 
have (no summation over i) 

A TE)y,^im^/D){ip\, i = 0,1, . . . , D - 1, D > 3. (5.14) 
Dkoiy - a) 

For D = 3 the radial stress {TI^)\j diverges logarithmically. In the case D = 2 the corresponding 
VEVs are finite on the boundary of the brane. 

Now we turn to the investigation of the general formulae for the brane-induced VEVs in the 
limit m/ku ^ 1- In this limit one has v ~ m/ko and the order of the modified Bessel functions 
in formulae (14. 5p and ()4.14p is large. Introducing a new integration variable x ux, we replace 
the modified Bessel functions by the leading terms of the corresponding uniform asymptotic 
expansions for large values of the order (see, for instance, [51]). By this way it can be seen that 
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the main contribution into the integral comes from the lower limit and to the leading order we 
find 



/ 2\ ^ k^-\z/za)^{m/kD)^l^~^ 1 - tanh(am) lu ^^ i , ^^ i...^ 

)^ " 2Dn^2^D,2^,2 /,2 _ i)D/2 1 + tanh(am) ["^("^Z^^) ^"(^/^")] ' ^^'l^) 

As we could expect, in this limit the VEV is exponentially suppressed. By the similar way, for 
the corresponding VEV of the energy-momentum tensor we obtain (no summation i) 

« (1 - 4e)m2(v,2)b, {TSh « ^(^0 )b. (5.16) 

The last relation is found most easily from equation (|4.17|) . 

In figures [T] and [2] we have plotted the dependence of the brane induced parts in the VEVs 
of the energy density and radial stress on z/ Za and ako for a minimally coupled massless scalar 
field = 0) in the case D = A. The first of these parameters is related to the distance from the 
boundary of the brane by the formula z/ Za = exp[A;£)(y — a)]. The curves for a = correspond 
to the VEVs in the RS 1-brane model. Recall that for a conformally coupled massless scalar 
field the brane induced VEVs vanish. Note that in D = 4 the conformal anomaly is absent and 
for massless scalar fields the VEV of the energy-momentum tensor in the free AdS spacetime is 
zero. 




Figure 1: The part in the VEV of the energy density, kj^ (^o)b5 induced by the brane as a 
function on z/ Za and ako for a minimally coupled massless scalar field in D = 4. 



5.2 Interior region 

Now let us consider the vacuum polarization effects inside the brane for the model with fiat 
interior. The corresponding eigenfunctions have the form given by Eq. (j2.8p with f\{y) = 
R{y, k, A) and the function R{y, k, A) is defined by formula (j5.4p . Substituting the eigenfunctions 
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Figure 2: The part in the VEV of the radial stress, fc^^ ^(r^)b, induced by the brane as a 
function on zj Za and akY) for a minimally coupled massless scalar field in D = 4. 



into the mode sum formula for the corresponding Wightman function one finds 



(0|(^(x)(^(x')|0) 

















dX 


^0 





dk 



(27r)^-i 

A cos{kyy) cos(fcyy')e"''^^'+^'^*"*'^ 



For the evaluation of the expression on the right we use the relation 

C{e^^^>'y<',Hl'\Xza)} 



(5.17) 



cos '^{kyo) 



irk 



D 



J^{XZa)+Y,^{XZa) 2ky 



^ s^,2 C{cos{kya), Hjf\Xza)} 

where 5s = (—1)*^^, and for given functions f{z) and g{z) we have introduced the notation 

'D 



(5.18) 



C{f{u),g{v)} = vf{u)g'{v) + 



u 



2 ^^^j^(^) ako 



(5.19) 



with u = kytt and v = Xza- 



By taking into account identity (I5.18|) . the Wightman function is written in the form 



{o\ip{x)ip{x'm 



{zakn)'^ / dk 



dX' 



(27r)^ 70 kyVk^T)? 



X cos{kyy) cos{kyy') 



1- E 



1 ^ C{e'^^''y'',Hl''{Xza)} 



2 ^^12 C{cos{kya), Hjf' (Xza)} 



(5.20) 
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To transform the A-integral with the first term in the square brackets, we present this integral 
as a sum of the integrals over the intervals (0, me~^°"') and (me~'^^", cxd). In the second integral 
changing the integration variable to ky we find 



kyy 



dX— cos(kyy) cos(kyy') = e / dky , 



X cos{kyy) cos{kyy') -i dX — —j===- cosh(y7/(A)) cosh(2/'r7(A)),(5.21) 

Jo r]{X)\/k^ + X-^ 



with the notation ?/(A) = V — A^e^'^o". In the first integral on the right of ()5.2ip . A is a 
function of ky given by formula (j5.4p . Let us consider the A-integral with the second term in the 
square brackets of ()5.20p . Note that in accordance with formula ()5.4p for ky, the corresponding 
integrand has branch point at A = me~^'^°'. For definiteness we will assume that this point is 
circled from above by a semicircle of small radius in the complex A-plane. With this choice, 
assuming that |y| + |y'| + |i — t'| < 2a, we rotate the integration contour in the complex plane 
A by the angle 7r/2 for s = 1 and by the angle —it/2 for s = 2. In the second case the integral 
with s = 2 is equal to the corresponding integral over the negative part of the imaginary axis 
plus the integral over the contour C depicted in figure El By taking into account that on the 
upper/lower parts of the contour C we have ky = zizirj{X), for this contribution one finds 

Xcos{kyy) cos{kyy')e-'^^^^+^^*-*''^ C{e-'''y'', Hl^\Xza)} 

UA ■ 



c 2kyVk^ + A2 C{cos{kya),Hi^\Xza)} 

dX cosh(y7?(A)) cosh(y'77(A)). (5.22) 

Now we see that the the integral over the contour C cancels the second term on the right of 
(I5.2ip . In the integrals over the imaginary axis the parts over the intervals {0,ik) and (0, —ik) 
cancel out. After introducing the modified Bessel functions, one can see that the Wightman 
function is presented in the form 



where 



{0\ip{x)^{x')\0) = Go{x,x') + Gi{x,x'), (5.23) 
Go{x, x') = {zakDf-^ J (^^^"'■^""" ^ dky ^yj^^ ^os{kyy) cosikyy'), (5.24) 



and 



(27r)^ J Jk C{cos\i{>i{x)a),Ku{xZa)] 

cosh(j<(x)y) cosh(x(x)y') 



{x)^/ x"^ — fc^ 



cosh[Vx2-fc2(t-t')]. (5.25) 



In (|5.25p . >c{x) = \J x^e^^oa _|_ ^^2^ g^j^^^j have used the notation (|5.19p with g(v) = Ky(y) and 
/(u) = e~",coshn for the numerator and denominator, respectively. 

In (j5.24p introducing a new integration variable K = ke'^^'^ we see that Go(x,x') is the 
Wightman function in the Minkowski spacetime in coordinates (j5.ip orbifolded along the y- 
direction: 

Gq{x,x')= I dK — r, — / dky , = cos{kyy)cos{kyy'), (5.26) 
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X 



c 



me 



-kr,a 



Figure 3: Integration contour C in the complex A-plane. 

where X = [X^^X"^,... ,X^^^) with coordinates X^ defined by relation (|5.2|) . This function 
differs by the factor 1/2 from the Wightman function for a plate in the Minkowski spacetime 
located at y = on which the field obeys the Neumann boundary condition. Evaluating the 
integrals, function (j5.26p is presented in the form 

1 , 

Go(x,x') = -Gm{x,x') + 



2 ™' ' ' 2(27r)(^+i)/2 

K(D-l)/2 (m V|X - X'|2 + (y + y'f - {X^ - X^'f 



(|X - X'|2 + (y + y'f - (XO - X0')2) 



2n(D-1)/4 



(5.27) 



where Gm{x,x') is the Wightman function for the Minkowski spacetime. Note that the expres- 
sion for Gm{x-, x') /2 is given by the same formula as the second term on the right of (|5.27p with 
the replacement y + y' y — y' ■ The term Gi (x, x') on the right of formula (I5.23P is induced by 
the AdS geometry in the region y > a. For a conformally coupled massless scalar field one has 
v = 1/2 and by using definition (I5.19P it can be explicitly checked that C{e~", K^{u/ak£))} = 0. 
Hence, in this case the part Gi{x,x') of the Wightman function vanishes. 

Now we turn to the evaluation of the renormalized VEV for the field square. The renor- 
malization corresponds to the omission of the part coming from the Minkowskian Wightman 
function in (j5.27p . As a result the VEV for the field square is presented in the form 

{v>')'^^ = {^')^:Sn + {^')t^- (5.28) 

Here the part (</'^)o reli is obtained from the second term on the right of (I5.27P in the coincidence 
limit and is given by the formula 



2y(mt) _ ^ ii:(g-i)/2 (2my) 

2(27r)(^+i)/2 (2my)(^-i)/2 



\'r /o,ren — ^^^_\( n^^^ /•) x('n-iW9 • {o.z,i)) 



The second term in the right hand-side of formula (|5.28p is obtained from (|5.25p taking the 
coincidence limit of the arguments and using the formula 

dk / dx , ^ ' = — — rr^Bl —,-] / dx x^+'~^g(x), (5.30) 

Jk V^^^ T{{D -l)/2) \ 2 '2jJo ^' ^ ^ 

with B{x,y) being the Euler beta function. This gives the result 

i^p't'^ = £ - cosh\xy)U^{x), (5.31) 
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with the notation 

U.{x) = (5.32) 
C{cosh(ax), K,y{\/x'^ — m? /ko)} 

This part in the VEV of the field square is induced by the AdS geometry in the exterior region. 
Note that if condition (15. Sp is satisfied the denominator in (15.32P is negative. Recall that under 
this condition there are no modes with imaginary A. 

The integral on the right of formula (j5.3ip is finite for < a and diverges on the boundary 
of the brane \y\ = a. In order to find the corresponding asymptotic behavior, we note that 
for points near the boundary the main contribution comes from large values x. By using the 
corresponding asymptotic formula for the function Ki,{x), to the leading order near the boundary 
y = a we find 

, 2amt) kpi^-^p) DT{{D -l)/2) 

W )b ~ -(4vr)(^+i)/2 (D-2)ia-y)D~2- (5.33) 

In the limit am ^ 1 the main contribution into the integral in formula ()5.3ip comes from the 
lower limit and one finds 



2\{mt) _ Bm^ ^cosh^(my) _2am r _ 1^/2 - 2^D + m/kp - . 



As we could expect, in this limit the VEVs are exponentially suppressed. For large values of the 
AdS curvature assuming that o, m l/ko, from definition (j5.19p it is easily seen that under 
the conditions ^ ^ l/(a/c£)), ^ — ^ l/(a^_D) the asymptotic behavior of ((/J^)!,™*^ is obtained 
substituting in formula ()5.3ip 

UAx) - ^2^- (5-35) 

The corresponding expression coincides with the VEV induced by Dirichlet boundary located 
at y = a in the Minkowski spacetime orbifolded along the y-direction with the fixed point 
y = 0. For special cases of minimally and conformally coupled scalar fields, in the limit under 
consideration the leading terms are obtained substituting in (j5.3ip 

f^-W « -^2^'^ = 0' (5-36) 
Uu{x) « -^-^=^^^ -, ^ = ^D. (5.37) 

In the minimally coupled case the corresponding limiting value coincides with the VEV induced 
by Neumann boundary located at y = a in the Minkowski spacetime orbifolded along the y- 
direction with the fixed point y = 0- 

Now let us consider the VEV for the energy-momentum tensor. As in the case of the field 
square, it is presented in the form 

{Tt)^^ = {ThSZ + {Tht"\ (5.38) 

where (Tifc)o is the vacuum energy-momentum tensor in the Minkowski spacetime in coordi- 
nates (|5.ip orbifolded along the y-direction and the presence of the part {Tik)^^^^ is related to 
that the geometry in the region y > a is AdS. For the first part one has (no summation over i) 

^^^^^"^'^ = (4^my)(^+i)/2 (2"^2/) (2e " 1) + (1 " 4^ myK^n+3)/2 (2my)} , (5.39) 
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with i = 0,1, . . . , D — 1 and (T'^)o = 0. For a massless scalar field this leads to the result 
(no summation over i) 

(int) _ D{C- gp) / D + l \ , . 

For the second term on the right of (j5.38p we find (no summation over i) 

-D/2 



^ cosh2(xy) + — 3^-^ 



cosh^(xy) - ^ 



(5.41) 



with i = 0, 1,...,-D — 1, and the function Uy{x) is defined by formula (|5.32p . Note that the radial 
stress inside the brane does not depend on spacetime point. This result could be also obtained 
directly from the continuity equation. As it has been mentioned before, for a conformally coupled 
massless scalar we have U^{x) = and, hence, the parts in the VEV of the energy-momentum 
tensor given by formulae (|5.4ip . (j5.42p vanish. From formula ()5.40p it follows that in this case 

the part (77)q°*^ vanishes as well. 

Now let us consider the VEV of the energy-momentum tensor near the brane core. In the 
way similar to that for the exterior region we find (no summation over i) 

(int) D-l ,^c,(int) DkD{C-(DfT{{D + l)/2) 

with i = 0,1, . . . , D — 1. For a conformally coupled scalar field the corresponding asymptotic 
behavior is given by formulae (j5.14p . In the limit am ^ 1 to the leading order one has (no 
summation over i) 

- ' /2 (4e-l)[2cosh^(my)-l], (5.44) 

,. ^, D D+l -2am 

2(47ram)-^/^ 

For large values of the AdS curvature, a,m ^> l/ko, the leading terms for the VEV of the 
energy-momentum tensor are obtained from formulae (j5.4ip . ()5.42p . by making substitutions 
(|5.35p for g » l/(aA;£)), g — g^i ^ l/(afc_D), and (j5.36p . (j5.37p for minimally and conformally 
coupled scalars. 

In figure H] we have plotted the dependence of the part in the VEV of the energy density 
induced by the exterior AdS geometry in the region inside the brane as a function of y/a and 
ako for a minimally coupled massless scalar field in the case D = 4. 

The corresponding radial stress does not depend on y and is depicted in figures [5l [6] as a func- 
tion of the parameters am and ako for minimally and conformally coupled fields, respectively. 
We recall that for a conformally coupled massless scalar field the corresponding VEVs vanish for 
both field square and the energy-momentum tensor. The perpendicular interior vacuum force 
acting per unit surface of the brane boundary is determined by —{TJ^)^^^\ As it is seen from 
figures O [6] for minimally and conformally coupled scalars these forces tend to decrease the brane 
thickness. 
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Figure 4: The energy density inside the brane, o (Tq)}^™ , induced by the AdS geometry in 
the exterior region as a function on ako and y/a for a minimally coupled massless scalar field 
in D = 4. 




Figure 5: The part in the radial stress inside the brane, a^^^ {TI^)\!^ , induced by the AdS 
geometry of the exterior region as a function on am and aki) for a minimally coupled scalar field 
in D = 4. 
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Figure 6: The same as in figure [5] for a confer mally coupled scalar field. 

6 Conclusion 

In braneworld models the investigation of quantum effects is of considerable phenomenological 
interest, both in particle physics and in cosmology. In the present paper we have considered 
the one-loop vacuum effects for a massive scalar field with general curvature coupling parameter 
induced by a Z2-symmetric thick brane on the {D + l)-dimensional AdS bulk. The previous 
papers on the investigation of the vacuum polarization by the gravitational field of the brane 
are mainly concerned with the idealized thin brane model, where the curvature has singularity 
at the location of the brane. Here we consider the general plane symmetric static model of 
the brane with finite thickness, described by the line element (|2.2p . Among the most important 
characteristics of the vacuum, which carry information about the internal structure of the brane, 
are the VEVs for the field square and the energy- momentum tensor. In order to obtain these 
expectation values we first construct the positive frequency Wightman function. In the region 
outside the brane this function is presented as a sum of two distinct contributions. The first one 
corresponds to the Wightman function in the free AdS geometry and the second one is induced by 
the brane. The latter is given by formula p.30p . where the tilted notation is defined by formula 
(|2.3ip with the coefficient from (j2.32p for the model without an infinitely thin shell on the 
brane boundary. This coefficient is determined by the radial part of the interior eigenfunctions 
and describes the infiuence of the core properties on the vacuum characteristics in the exterior 
region. In the case of the model with a thin shell on the boundary of the brane, the derivatives 
of the metric tensor components are discontinuous on the brane surface. This leads to the 
delta function type contribution in the Ricci scalar and, hence in the equation for the radial 
eigenfunctions in the case of the non-minimally coupled scalar field. As a result, the radial 
eigenfunctions have a discontinuity in their slope at the brane boundary. This leads to an 
additional term in the coefficient of the tilted notation which is proportional to the trace of the 
surface energy- momentum tensor (see Eq. (jS.lOp ). 

By using the formula for the Wightman function, in section [3] we have investigated the 
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influence of the non-trivial internal structure of the brane on the VEVs of the fleld square 
and the energy- momentum tensor. The parts in these VEVs induced by the brane are directly 
obtained from the corresponding part of the Wightman function for the case of the fleld square 
and by applying on this function a certain second-order differential operator and taking the 
coincidence limit for the energy-momentum tensor. In the general plane symmetric model for 
the brane interior these parts are given by formulae (j4.2p and (j4.10p for the fleld square and 
the energy-momentum tensor respectively. These formulae are further simplifled for models 
with Poincare invariance along the directions parallel to the brane taking the form (j4.5p . (j4.14p . 
For a conformally coupled massless scalar field the corresponding energy-momentum tensor 
vanishes. The parts in the VEVs of the field square and energy-momentum tensor induced by 
the brane diverge on the boundary of the brane. The surface divergences in the VEVs of the 
local observables are well-known in quantum field theory with boundaries and are investigated 
for various boundary geometries. At large distances from the brane the brane induced VEVs 
are suppressed by the factor e~'^^^^^ . In the limit of strong gravitational fields corresponding 
to large values of the AdS energy scale ko-, for points not too close to the brane the parts in the 
VEVs induced by the brane behave as j^^^^^-'^'^'^Diy-a) with upper/lower sign corresponding to 
the energy- momentum tensor /field square. In this case the relative contribution of the brane 
induced effects are exponentially suppressed with respect to the free AdS part. 

As an application of the general results, in section [5] we have considered a simple model with 
fiat spacetime in the region inside the brane. The corresponding surface energy-momentum 
tensor on the boundary of the brane is obtained from the matching conditions and has the 
form given by Eq. (15. 3p . The brane induced parts of the exterior VEVs in this model are 
obtained from the general results by taking the function in the coefficient of the tilted notation 
from Eq. (j5.7p . We have also investigated the vacuum densities inside the brane. Though the 
spacetime geometry inside the brane is Monkowskian, the AdS geometry of the exterior region 
induces vacuum polarization effects in this region as well. In order to find the corresponding 
renormalized VEVs of the field square and the energy-momentum tensor we have presented the 
Wightman function in the interior region in decomposed form (j5.23p . In this representation the 
first term on the right is the Wightman function in the Minkowski spacetime orbifolded along 
the direction perpendicular to the brane and the second one is induced by the AdS geometry in 
the exterior region. The corresponding parts in the VEVs of the field square and the energy- 
momentum tensor are given by formulae (j5.3ip . (j5.4ip . (j5.42p and are investigated in various 
asymptotic limits of the parameters. For a massless conformally coupled scalar field these parts 
vanish. In the general case of the curvature coupling parameter, the corresponding radial stress 
is uniform inside the brane and determine the interior vacuum forces acting on the boundary of 
the brane. The results of the numerical calculations plotted in figures [5l [6] show that for both 
minimally and conformally coupled scalar fields these forces tend to decrease the thickness of 
the brane. When the brane thickness tends to zero, from the formulae of the model with fiat 
interior the corresponding results in the RS 1-brane model are obtained. 
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